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Abstract 

The calculation of the quartic mass terms at order a^ to the hadronic R ratio is 
described. It is based on the operator product expansion for the quark current 
correlator, combined with an application of the rcnormalization group equation. 
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Abstract 

The calculation of the quartic mass terms at order a^ to the hadronic R ratio is 
described. It is based on the operator product expansion for the quark current 
correlator, combined with an apphcation of the renornialization group equation. 



1. Introduction 

The hadronic R ratio is one of the most important 
quantities of QCD. It is defined as the total cross 
section for hadronproduction in e~^e~ annihilation, 
compared to the production rate for muon pairs: 
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Its inclusive character is a welcome feature both 
from the experimental and the theoretical point of 
view. For example, the calculation of R{s) can be 
performed fully perturbatively, as long as bound 
state effects which emerge close to threshold can 
be neglected. 

The standard way to compute the hadronic R 
ratio is to apply the optical theorem which relates it 
to the imaginary part of the polarization function: 



R{s) = 127rlmn(s + ie) . 



(2) 



Examples for diagrams contributing to 11 (g^) 
up to order a^ are shown in Fig. m. The 
analytic evaluation of such diagrams is only 
possible up to two-loop order, if a non-vanishing 
quark mass is assumed. At three-loop order, 
only the massless case can be computed fully 
analytically. Nevertheless, mass effects at 0{a1) 
could be evaluated to very high accuracy, using 
two different approximation methods. One is based 
on asymptotic expansions ||l|, yielding the result in 
terms of a series in m? jq^ times logarithmic terms 
of the form ln(— q^//i^) and ln(— q^/rn^), where 
q^ is the renormalization scale 0. The second 
method uses the moments of n(q^), defined as 
C„ = (a"/a(g2)")n(g2)|^2^o, to reconstruct I{{q^) 
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Figure 1. Sample diagrams contributing to n(q^) up 
to four-loop level. 



in the full complex plane. The result for R{s) is 
therefore valid over the full energy range in this case 

i- 

At ©(a^), however, either of the two methods 
above inevitably leads to four-loop diagrams whose 
evaluation is possible only for a limited number 
of special cases at the moment. Nevertheless, 
besides the massless approximation [Q which could 
be reduced to a three-loop calculation using a 
method called infrared re-arrangement [pi , also the 
uadratic mass corrections at 0{a'l) are available 

They were obtained by a pure renormalization 
group (RG) analysis using the terms at 0{a1) as 
input information. 

At lower order in as, the quartic {m^) mass 
corrections turn out to be quite sizable and are 
thus needed to provide a good approximation to 
the full result over a reasonably large energy range. 
This was the main motivation for us to evaluate 
these corrections at 0{a'l). The method which was 
followed has been developed in Q. At that time it 
was used to determine for the first time the quartic 
corrections at 0{a1). The strategy is to evaluate 
the coefficient functions and the matrix elements of 



the operator product expansion (OPE) for n(g^) to 
(in our case) the three-loop approximation and to 
use the RG equation in order to extend the vahdity 
of the result to 0{a^). 

2. Operator product expansion (OPE) 

The polarization function is defined as the vacuum 
expectation value of the time ordered product of 
two currents, j^ — V'T/jV'j where '0 be a quark field 
with mass m. In the limit of large q^, n(g^) can be 
written asymptotically as an OPE: 
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Figure 2. 
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Sample diagrams contributing to (d) and 







where (On) = (OlOnlO). Only Lorentz scalar 
operators contribute due to the external vacuum 
states. The function n(q^) in (g) is defined as the 
transversal part of 11^^ (q) , 



n(g2) 



1 1 



-ffM" 



9 



n^.Cg)- (4) 



An analogous definition is used to define transversal 
coefficient functions C„. The operators up to 
dimension four are given by 
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(the superscript "(4)" will be dropped in what 
follows). In addition to the renormalization of 
m and a^, one has to take into account operator 
mixing according to 
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In the MS scheme which we will use throughout 
this paper, operators of different mass dimension 
do not mix under renormalization. For the 
dimension-4 operators the renormalization matrix 
Znm has been evaluated in g]. It was expressed 
through the renormalization constants for m and 
Us, plus the one for the QCD vacuum energy, 
Zq. If no normal ordering is applied to the 
operators, the perturbative expressions for their 
vacuum expectation values are not necessarily zero. 
They lead to tadpole diagrams, examples of which 
are shown in Fig. 0. The matrix elements for the 
"trivial" operators are given by (m^") = m?''^ to all 
orders of perturbation theory. The ones for Oi and 
O2 can be obtained up to 0{a1) by performing a 



Figure 3. 

and C3. 



Sample diagrams contributing to Ci, C2, 



three-loop calculation, as is evident from a look at 
the diagrams in Fig. g. 

Concerning the coefficient functions it has 
been shown in P] that they depend only on 
the external momentum if minimal subtraction is 
used as the renormalization scheme. The most 
elaborate method to reduce their evaluation to 
massless propagator diagrams is called the "method 
of projectors" iQ. It leads to diagrams like 
the ones shown in Fig. 0, where it is understood 
that certain derivatives w.r.t. the quark mass and 
the momentum carried by the external quark and 
gluon lines are to be applied to the integrand. 
Before integration, however, all masses and external 
momenta (except q) can be set to zero. A glance at 
these diagrams shows that a three-loop calculation 
yields Ci and C2 with ©(a^), but C3 only with 
0{al) accuracy. Furthermore, one realizes that the 
imaginary parts arising from Ci and C2 start to be 
non-zero only at 0{as)- Thus the a^ calculation for 
the matrix elements mentioned before is sufficient 
to determine the imaginary part of the products 
Gi{Oi) andC2(02) toO(a3). 

The crucial step to reconstruct the a^ terms 
of C3 is now to employ RG invariance of the 
dimension-4 piece of the OPE: 
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The effect of A^^g4^ on the On is determined by (||): 



H,,2^" ~ 7 ,lnmOn 



dfj, 



(8) 




Figure 4. Quartic mass corrections to charm pair 
production. 



where 7mn is the anomalous dimension of the 
operators which is related to the renormalization 
matrix Z„„j in the standard way. For the coefRcicnt 
functions, the derivative w.r.t. /i will be rewritten 
as 
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with the QCD (3 function and L = ln(— /i^/g^). 
Solving (0) for ^Cs and keeping only terms 
proportional to O3, one obtains 



^C3 = -a./3-^C3-47™C3 
oL oas 
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where 733 = 47™ has been used (7™ is the quark 
anomalous dimension). Since /3 and 7™ do not carry 
terms of 0{a'^), the r.h.s. of (M) is known to 0{al), 
and so is the l.h.s. The logarithmic terms (and thus 
the imaginary part) of C3 may therefore be obtained 
by trivial integration. 

3. Results 

Let us consider the case of charm quark pair 
production. The general result for the terms of 
order a'^m^/s'^ to R{s) will be given elsewhere 
Jll| . The upper limit of the energy region we are 
interested in is dictated by the bb threshold. The 
lower one is such that, on the one hand, the massless 
charm approximation fails, and on the other hand, 
one is sufficiently far above the cc threshold in order 
for the small mass expansion to be valid (see Ql for 
a detailed analysis on the convergence properties at 
0{al)). Fig. ^ shows the result for the contribution 
where the external photon couples directly to the 
charm quark. There is also a contribution arising 
from diagrams with only massless quarks coupling 



to the photon and the charm quarks emerging only 
through their coupling to gluons. In turns out, 
however, that the latter piece is tiny, so we will 
neglect it here. The notation is given by the 
following decomposition: 



Rc{s) = 3 rc,o(s) + rc,2is) + r^iis) + 



(11) 



which means that r'c.o(s) corresponds to the 
massless approximation and rc,2(s) and rc,4(s) 



contains the mass corrections of order 



and 



respectively. 

The charm quark mass is taken in the MS 
scheme and the renormalization scale is set to 
fi^ = s. One observes that the higher order 
terms in as are all of comparable magnitude, and 
the convergence of the perturbative series is not 
manifest. Compared to the Born result, however, 
the corrections in a^ are small such that the 
stability of the full prediction is almost unaffected 
by the higher order terms. 
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